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I. INTRODUCTION 

Understanding the phase diagram of strongly inter- 
acting theories will unveil a large number of theories of 
fundamental interactions useful to describe electroweak 
symmetry breaking, dark matter and even inflation [T- 
HJ. To gain a coherent understanding of strong dynamics 
besides the SU(N) gauge groups 0161, one should also 
investigate the orthogonal, symplectic and exceptional 
groups. SO(N) and SP(2N) phase diagrams were investi- 
gated with analytic methods in 0, while the exceptional 
ones together with orthogonal gauge groups featuring 
spinorial matter representations were studied in J5J. So 
far lattice simulations have been mostly employed to ex- 
plore the phase diagram of SU(N) gauge theories while 
a systematic lattice analysis of the smallest symplectic 
group was launched in 1191 . 

Here we move forward by analyzing on the lattice 
the dynamics of the SO(4) gauge group with two Dirac 
fermions in the vector representation of the group. This 
choice is based on the following theoretical and phe- 
nomenological considerations. The theory is expected 
to be below or near the lower boundary of the conformal 
window [TO] , and therefore break chiral symmetry. 
The theory can be used as a technicolor IfTTl [12 1 tem- 
plate similar, from the global symmetry point of view, to 
Minimal Walking Technicolor (MWT) I51H31I141. How- 
ever, the fermion representation is such that, differently 
from MWT, one cannot construct composite fermions 
out of one techniquark and one techniglue. This re- 
moves immediately the presence of fractionally charged 
states with the simplest choice of the hypercharge as- 
signment. Furthermore the technicolor theory leads 
to a weak isotriplet with the neutral member being an 
ideal dark matter candidate JT) [lOl , the isotriplet Tech- 
nicolor Interactive Massive Particle (iTIMP). This state 
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is a pseudo Goldstone and therefore can be light with 
respect of the electroweak scale making it an natural 
candidate to resolve some of the current experimental 
puzzles IflOl |l5l . The first model featuring composite 
dark matter pions appeared in 1 16 17] and the first study 
of technipion dark matter on a lattice appeared in [9J. 

Due to the reality of the fermion representation the 
quantum global symmetry group is SU(4) expected to 
break spontaneously to SO(4), yielding nine Goldstone 
bosons. Once gauged under the electroweak theory 
three are eaten by the SM gauge bosons. Six additional 
Goldstone bosons form an electroweak complex triplet 
of technibaryon with the neutral isospin zero component 
to be identified with the iTIMP of |TOl. 

SO(4) is a semi simple group, SO(4) = SU(2)®SO(3), 
and it has a non-trivial center Z2. The theory is asymp- 
totically free and since the two-loop |3-function for dif- 
ferent number of flavors looses the infrared zero for 
Nf = 2.3 while the all-orders beta function 11811191 pre- 
dicts the anomalous dimension of the mass to be unity 
for Nf = 2.86 we expect that chiral symmetry breaks for 
two Dirac flavors. However, we want to confirm here 
this result via first principle lattice simulations. Further- 
more there is also the possibility that the theory shows a 
certain degree of walking |[7l [T0ll20H22l unless the phase 
transition is of jumping type |23, 24 1. Jumping confor- 
mal phase transitions have been demonstrated to occur 
in a wide class of theories 1 25 1 . 

As a natural first step, we study the phase diagram in 
the (jS, ffio)-plane to find the relevant region of parame- 
ter space to simulate. We then determin the zero PCAC 
mass line as well as the strong coupling bulk phase tran- 
sition line. In addition, we report on the pseudoscalar, 
vector and axial vector meson masses. From the mea- 
sured spectrum we infer that the theory breaks chiral 
symmetry dynamically. Part of these results appeared 
in (261. 

In Section [IT] we present the analytic expectations for 
the phase diagram of SO(N) as function of the number 
of flavors. We also summarize the expected breaking 
pattern of the quantum global symmetries for theories 
below the conformal window. We also prove the spec- 
tral degeneracy between certain diquarks and ordinary 



meson-like states. In Section [ITT] we recall the lattice for- 
mulation of the theory and summarize the physical ob- 
servable studied here. The results of the simulations are 
reported in Sect. [TV| and conclude in Sect. |V| 



II. ORTHOGONAL CONFORMAL WINDOW AND 
CHIRAL SYMMETRY BREAKING PATTERN 

The two loop |3-function for an SO(JV) theory with Nf 
Dirac f ermions transforming according to the vector rep- 
resentation of the gauge group is 
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A naive estimate of the lower bound of conformal win- 
dow is given when the second coefficient b\ changes 
sign. For SO(4) this happens when Nf = || - 2.3. The 

corresponding values for three and four-loops in the MS- 
scheme are Nf — 1.8 and Nf — 3.0. The all-orders beta 
function predicts as lower boundary Nf = 2.86, see Fig. [TJ 
Hence, perturbative and nonperturbative methods sug- 
gest that chiral symmetry breaks for two Dirac flavors. 
However, lattice simulations can seal this expectation. 
Since the vector representations of orthogonal groups are 
real the quantum global symmetry of the theory is, for a 
generic Nf SU(2N f ) which is larger than SU(Nf)®SU(N f )® 
UV(1) valid for complex fermion representations. The 
reality property of the representation translates in the 
following property of the Dirac operator 



(U> + m)Cy 5 = Cy 5 (U> + m)\ 



(3) 



where ty = y''(<? f , - igA a „T a ), a = 1, ...,d[G] where d[G] is 
the dimension of the adjoint representation of the gauge 
group and C = iy a y 2 is the charge conjugation operator. 

The global SU(2Nf) is assumed to break to the maximal 
diagonal subgroup 



SU(2N" f ) 



SO(2N f ) , 



(4) 



for the massless theory and for Nf below the conformal 
window. A common mass for the Dirac fermions leads 
to the same pattern of explicit symmetry breaking. The 
explicit interpolating operators for the Goldstones can 
be naturally divided in three independent antifermion- 
fermion bilinears 



{\, f y ip f , , 



with f and /' the flavor indices / = 1, 2 and six difermion 
operators 



tfCy^f, 



and 4>fy 5 C4>J, . (6) 
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FIG. 1: Conformal window of SO(N c ) with N f Dirac fermions 
in the fundamental representation. Upper bound is when 
asymptotic freedom is lost. Lower bounds are 2-loop (red, 
dashed), 3-loop (yellow, dotted) and 4-loop estimates (green, 
solid). 



The reader can find a useful summary of the global sym- 
metry breaking patterns tailored for lattice computations 
in [27 1 while applications to beyond standard model 
physics for similar patterns appeared in |16 28]. No- 
tice that whereas the usual pions have odd parity, the 
corresponding diquarks are parity even. It was noticed 
in that when fermions are in a pseudoreal represen- 
tation, the diquark correlator is exactly identical to the 
corresponding mesonic correlator. In appendix [A] we 
give a similar proof for fermions in real representations. 
The proof uses the symmetry |3} of the Dirac operator 
along with the y 5 -hermiticity y^(W + m)y 5 = + mf 
property. The result can be stated as 



(T) 

where cy^ is the correlator for the operator ipfTijjf, 
and cj^ t is the correlator for the corresponding di- 
quark operator ip T f TCip p . T can be any of the matrices 



,(*-y) 



cv.,,(x-y), (7) 



(5) l,y 5 ,yf',yf'y 5 . 



Having discussed the generic features expected for 
orthogonal groups we now turn to the lattice formulation 
and results for the relevant case of SO(4) with two Dirac 
flavors. 
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III. LATTICE FORMULATION 

In this work we have used the Wilson prescription for 
the lattice action 



S = S F + S G , 



where 



X fi,V<ff c 



(8) 



(9) 



is the Yang-Mills gauge action. We have normalized the 
lattice spacing to a - 1. U„ v (x) is the plaquette defined 
in terms of the link variables as 

U flv (x) = U lt (x)U v (x + p)U T Jx + fi + v)Ul(x + v) . (10) 
The Wilson fermion action is 

(11) 



f x -y 



with f running over fermion flavors and the Wilson- 
Dirac matrix M(x, y) given by 

Y u M(x,y)^(y) = (4 + m )xjj(x) 
y 

- \ E [ (1 + - aw* ~ w (12) 

+ (l-yf l )U (l (x)il>(x + £L)]. 

Here the gauge and spinor indices have been suppressed. 
The bare parameters are the inverse of the bare coupling 
jS = 2N c /gy appearing in the gauge action and the bare 
mass nto of the Wilson fermions. 

We employ the Partial Conservation of the Axial Cur- 
rent (PCAC) relation to define the physical quark mass 



mpcAC = lim 



1 d t V FS 



(13) 



t— >oo 2 X^pp 
where the currents are 

Vps(xo) = a 3 Yj (M^yoysM^ii^rsM )) • 

,Tl,l'2,.V3 

Vpp(x ) = a 3 Y (MxyysMxWiMysM )) ■ ( 14 ) 



Xl,X 2 ,X 3 



The meson masses are estimated using time slice av- 
eraged zero momentum correlators 

Cf if2 (x ) = a 3 Y K([HxW2(x)T4>iWM0)) . 

X-i,X2,X 3 

(15) 

where T = y 5 for pseudoscalar, T = (k = 1,2,3) for 
vector, and ysyk for axial vector meson. 
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FIG. 2: Lattice phase structure outlined on an 8 3 x 16 lattice. 
Circles represent points of critical bare mass where ?«pcac = 0. 
The transition between the bulk phase is of first order. The error 
bars represent the interval over which the measured average 
plaquette jumps. 



Volume 



Iterations Thermalization 



4.1,4.2. . . 4.9, 5.2, 5.4, 5.6 2000 500 

4 4.5 5, 5.5, 6, 7 5000 2000 

5.5, 7 5000 1500 

7 850 - 2000 600 



8 3 x 16 

12 3 x 64 
24 3 x 64 



TABLE I: Simulation parameters and thermalization times. For 
each coupling we performed multiple simulations with appro- 
priate bare masses. The thermalization column refers to the 
number of discarded initial configurations. 



IV. RESULTS 

The simulations were performed on three different lat- 
tices 8 3 x 16, 12 3 x 64 and 24 3 x 64 where in all cases the 
larger dimension is the temporal one. All the simula- 
tions were started from a random configuration and the 
first 500-2000 iterations were discarded. This is enough 
to thermalize the system for the quantities we measured. 
For a complete list of the simulations see Table [T] where 
we have omitted the values of the bare masses. 

The smallest lattice was used for exploration of the 
parameter space spanned by the bare mass m§ and the 
coupling jS. Fig. |2]shows an outline of the lattice phase 
structure measured on this 8 3 x 16 lattice. For small val- 
ues of jS the system is in a bulk phase not connected to 
continuum physics. The bulk phase is separated from 
the small coupling (large /3) phase by a first order phase 
transition. Fig. [3] shows the discontinuous behavior of 
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FIG. 3: Average plaquette (P) vs. /3 on an 8 3 x 16 lattice at three 
different values of the bare mass. 



FIG. 4: /Wpcac in units of inverse lattice spacing at three different 
couplings. The measurements are performed on a 8 3 x 16. 



the average plaquette when crossing the bulk phase tran- 
sition, for three different values of ffJo . The uncertainty on 
the location of the bulk phase transition shown in Fig. [2] 
is due to taking discrete values of B between simulation 
points. 

We can compare our result for the location of the bulk 
transition to previous studies of SO(N) pure gauge the- 
ories. Earlier simualtions focused mainly on the SO(3) 
gauge group E9l with the exception of [30J where also 
other values of N were considered. For SO(4) the au- 
thors of 1 30 1 find that the bulk phase transition happens 
for 4.62(3) < B < 4.87(3), which is in agreement with our 
result in Fig. [3] 

The critical line where the physical quark mass van- 
ishes is determined from the PC AC relation (13) . The 
critical line of rtiq = in the phase diagram (Fig.j2Fis con- 
structed by linear fits to the PCAC mass. Fig/Bl shows 
the bare mass dependency of the PCAC mass at three 
different couplings on the 8 3 x 16 lattice. 

A. Finite size effects 

According to the perturbative estimates discussed in 
section |ll] the running of the gauge coupling is expected 
to be slow. This also suggests that attention should be 
paid to finite size effects, which need to be estimated 
non-perturbatively by measuring physical observables 
as a function of lattice size. 

In the case of SO(N) pure gauge theories [29j |30l the 
bulk phase transition occurs at such a weak coupling that 
extremely large lattices are required for simulations in 



the confined phase, the one connected to the continuum 
physics. However, in the presence of dynamical quarks, 
we find that somewhat smaller volumes (24 3 X 64) are 
enough to probe the chiral regime of the system. 

In figure Fig.|5]the mass of the pseudoscalar meson nips 
and the PCAC quark mass wzpcac is plotted for different 
lattice sizes. The PCAC mass has little dependence on 
the lattice size being a UV quantity. The pseudoscalar 
meson mass, on the contrary, is very sensitive to finite 
size effects even if it is still somewhat heavy at the bare 
mass used in Fig. [5] 

Another interesting property which occurs in this 
model is a novel phase separation at small volumes. 

We observe the coexistence of two distinct phases 
which can be characterized by the average value of 
Polyakov loops wrapping around the three spatial di- 
rections taken on each time slice separately. In detail, 
the operators we consider are defined as 

L ^={m^k Tvl ^ m ' x) )' (16) 

where i + j + k are spatial directions. Fig. [6] shows the 
time resolved Polyakov loops on a 12 3 x 64 lattice at B = 7 
and nio = -0.3. The averaging is performed over 700 
configurations belonging to the same simulation. The 
coexistence of two phases with different values of L2 is 
clear from the figure. The phenomenon appears in all 
simulations performed on small lattices. The location 
of the phase boundaries and the direction in which the 
Polyakov loop has non-zero average is random. In some 
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FIG. 5: Finite size effects on m P s and tmpcac in units of inverse 
lattice spacing. The measurements are performed on a 24 3 x 64 
lattice at $ = 7 and m = -0.2. 
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FIG. 6: Average Polyakov loops wound around the three spa- 
tial dimensions computed at each timeslice of the lattice. This 
measurement was performed on a 12 3 x 64 lattice at p = 7 and 
m = -0.3. The values are averages over 700 configurations 
starting at 1800 where the system does not appear to thermal- 
ize further. 



cases more than two phase boundaries appear in the 
same system. 

The coexistence of two phases with different mesonic 
correlation lengths separated by a domain walls could 
explain the unusual decorrelation of mesonic operators 
observed for volumes smaller than 24 3 x 64. This is 
reflected in the rise of the effective mass plateaux of the 
pseudoscalar meson shown in Fig. [7] 

In order to understand whether these phase separa- 
tions are related to the presence of dynamical fermions 
we have also performed pure gauge simulations on 
12 3 X 64 lattices. The phase separation occurs also for 
the pure gauge. Thus the phenomenon seems to be a 
feature stemming from the pure gauge sector. 

We will not explore this feature further, but it would 
be interesting to continue its investigation in the future. 

In order to avoid the complications stemming from 
the phase separation described above we use 24 3 x 64 
lattices for the rest of the paper. 



B. Spectrum and chiral symmetry breaking 

We address the dynamical fate of the chiral symme- 
tries of the theory by determining the pseudoscalars and 
(axial) vectors spectrum. 

Fig.[8]shows the pseudoscalar, vector, and axial vector 
meson masses measured on a 24 3 X 64 lattice at f> = 7 
as the bare quark mass is decreased towards the critical 
value. At the lightest quark mass the pseudoscalar me- 



son has a mass of about mps — 0.15 in lattice units. This 
means that mpg • L 3.6 which is where finite volume ef- 
fects starts to become relevant. At large quark masses the 
vector and pseudoscalar are degenerate with the com- 
mon mass increasing linearly with the quark mass. At 
smaller masses the vector meson becomes heavier than 
the pseudosclar. This is consistent with dynamical gen- 
eration of a chiral scale. To see this more clearly the ratio 
of the vector and the pseudoscalar masses have been 
plotted in Fig. [9] Indeed the mass ratio approaches unity 
for large quark masses. However, when approaching the 
chiral limit the ratio increases signaling chiral symmetry 
breaking. In fact this result is consistent with the expec- 
tation that if spontaneous symmetry breaking occurs the 
vector meson remains massive whereas the pseudoscalar 
meson is massless. A diverging ratio frzy/ffips therefore 
indicates chiral symmetry breaking. This is the trend we 
observe in Fig. [9] However to nail this conclusion more 
studies have to be performed. 

The axial mass in the chiral limit is poorly determined 
Fig. [8] In the future we plan on improving its determi- 
nation. We will then be able to use it to infer interesting 
properties of the chiral transition. For example one can 
investigate whether the axial remains (near) degenerate 
with the vector in the chiral regime which could signify 
the onset of walking dynamics |28 3T[. 

To extract further properties of the theory we analyze 
in more detail the functional dependence of the pseu- 
docalar mass on the quark mass. It is well known that, 
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FIG. 7: Effective mass of pseudoscalar meson for two different 
volumes. 



FIG. 8: Pseudoscalar, vector, and axial vector meson masses 
measured on a 24 3 x 64 lattice at j8 = 7. 



meson fit 

ps chiral 

ps conformal 

ps alt. 1 
ps alt. 2 

vector chiral 

vector conformal 
vector alt. 1 

vector alt. 2 



fit function best parameter ^p/dof 



a ^Jm 

am 

a + bm 
a + byfin 

a + bm 

am 
a -\fm 

a + byfin 



a = 1.167(6) 
a = 4.69(3) 
a = 0.111(6) 
b = 2.9(1) 
a = -0.001(10) 
b = 1.17(4) 
a = 0.16(1) 
b = 2.3(2) 
a = 4.91(3) 
a = 1.231(6) 
a = 0.07(2) 
b = 0.96(7) 



0.43/2 
364/2 

6.4/1 
0.41/1 

3.3/1 

273/2 
18/2 

0.69/1 



TABLE II: Different types of fit functions in the chiral regime 
for the data with m identified with the m PC AC- 



for this kind of theories, spontaneously broken chiral 
symmetry leads to the Gell-Mann-Oakes-Renner rela- 
tion 1321 



m 



PS 



Km 



PCAC / 



(17) 



valid in the chiral limit, where A = -2(i/>i/>)//p g is a dy- 
namically generated scale. For conformal theories the 
behavior is different 15511341 . In [34 1 it was also shown 
that the instanton contributions to conformal chiral dy- 
namics can be neglected when the anomalous dimension 
of the mass operator is less than one. This property has 
been investigated and confirmed via lattice simulations 
in 1 35 J . A clever separation of the ultraviolet and infrared 
modes presented in 1 36 37 ] led to a better understanding 
of the conformal chiral scenario but without discussing 



the instanton contributions [34 J. Building upon these 
results an interesting method to determine the anoma- 
lous dimension of the fermion masses was put forward 
in 1581 . To sum up, for a conformal scenario the dy- 
namical scale A mutates into a fermion-mass dependent 
quantity 1541 and therefore 



In Fig. 



10 



must vanish as ?Wp CAt 
we plot the ratio wz ps /mpcAC ror decreasing 
fermion~mass. We see that the ratio approaches a con- 
stant for vanishing fermion masses which is consistent 
with the chiral symmetry breaking scenario $F7\ . 

In table In] we report the fit to the data for the depen- 
dence of the pseudoscalar mass as well as the vector 
mass as function of the ffipcAC within the believed chiral 
regime of the theory. This corresponds to the three low- 
est values of /«pcac where the ratio ni'L/mpcAC becomes 



roughly constant as shown in Fig. [10] The data points 
used for the chiral fits in the table are shown in Fig. 11 
The best fit curve, determined by the lowest ^ 2 /dof, for 
the pseudoscalar mass corresponds to the first line of the 
table which is in agreement with the GMOR expectation. 
It is remarkable that by even allowing for an offset of the 
mass value in the chiral limit the best fit demands the 
offset to vanish, see the last line of the table. We have 
tried also to test the possibility that the pseudoscalar 
mass vanishes linearly with the fermion mass and found 
that this is highly disfavored. If the theory would have 
been conformal we would have expected this case to fit 
much better. 

Similarly, by fitting the vector masses dependence on 
the fermion mass, in the lower part of table |n] we ob- 
serve a reasonable agreement with the expected chiral 
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FIG. 9: Ratio between pseudoscalar and vector meson masses 
measured on a 24 3 x 64 lattice at p = 7. 
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FIG. 11: The chiral fits to the pseudo scalar and vector meson 
masses on a 24 3 x 64 lattice at fi = 7. 




FIG. 10: Psudoscalar mass squared divided by the quark mass 
measured on a 24 3 x 64 lattice at § = 7. 



behavior of the theory. The two best fits correspond to 
the first and last line of the lower part of the table. We 
would have expected the first line to yield a better fit if 
chiral symmetry breaks like in ordinary Quantum chro- 
modynamics. We believe that for this case more statistics 
is needed to resolve which of the two cases is actually 



realized given that the data cannot yet differentiate be- 
tween the two. As for the pseudoscalar case the would 
be conformal case is highly disfavored (see second line 
of the lower part of the table). 

Using the identity for the hadronic correlators |7) we 
can immediately infer the baryonic diquark masses. 



V. CONCLUSIONS 

Orthogonal lattice gauge theories with dynamical 
fermions have so far been terra incognita. However, as 
explained in the introduction, these theories can be rel- 
evant for models of dynamical electroweak symmetry 
breaking as well as for the construction of interesting 
dark matter candidates. Furthermore to have a deeper 
understanding of strong dynamics it is essential to gain 
information on different gauge theories. We have cho- 
sen to start investigating the orthogonal gauge groups 
dynamics with a phenomenologically relevant example, 
i.e. the SO(4) gauge theory with two Dirac flavors trans- 
forming according to the vector representation of the 
group. 

We have uncovered the lattice phase diagram and 
shown that there is a novel phase separation phe- 
nomenon at small volumes which persists even in the 
pure gauge case. We have shown that the phase sepa- 
ration can be circumvented and the chiral regime of the 
theory studied using large but still feasible lattices. 

Finally we investigated the spectrum of the theory for 
the pseudoscalar, vectors and axial vectors. The results 
for the spectrum are consistent with chiral symmetry 



8 



breaking and strongly disfavor a conformal behavior. 

Appendix A: Diquark correlators 

A generic mesonic correlator will have the form 
c£>,(* - y) = Tr ([4>(x)Tf (x)f $(y)Tf(y)) , (Al) 
and the baryonic diquark correlator will have the form 

cJJ,(x - y) = Tr ([</> T \x)CTf (x)f f \y)CTf (y)) . (A2) 

Rewriting the diquark correlator slightly gives 

c^,(x -y) = Tv(T^(y)4>'(x)y T f C f (Y°) T [i/Wtof c) . 

(A3) 

Now we can invoke two identities 

(yff = -Cy^C 4 , (A4) 
V(*)ip(y) = C + ^(y)^(*)] r C. (A5) 



The latter identity follows from the symmetry of the 
Dirac matrix given in |3) along with y 5 -hermiticity 
y 5 (# + m)y 5 = (ty + mf . The identity jA5) extend to the 
Wilson lattice formulation of the Dirac matrix. This is 
demonstrated for pseudoreal representations in the ap- 
pendix of 0. Invoking the identities in the expression 
for the diquark correlator ( | A3[ ) we have 



cj>,(x - y) = Tr (r^'(y)^'(x) ? /°rV>(x)^(y)) = cJJ,(* - y) . 

(A6) 



A similar derivation holds for the antiparticles leading 
to the identity 0. 
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